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Abstract 

We investigate the scaling properties of the spin interfaces in the Ashkin- Teller model. These 
interfaces are a very simple instance of lattice curves coexisting with a fluctuating degree of 
freedom, which renders the analytical determination of their exponents very difficult. One of 
our main findings is the construction of boundary conditions which ensure that the interface still 
satisfies the Markov property in this case. Then, using a novel technique based on the transfer 
matrix, we compute numerically the left-passage probability, and our results confirm that the 
spin interface is described by an SLE in the scaling limit. Moreover, at a particular point of the 
critical line, we describe a mapping of Ashkin- Teller model onto an integrable 19-vertex model, 
which, in turn, relates to an integrable dilute Brauer model. 
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1 Introduction 

The scaling limit of interfaces in statistical models has attracted much attention since the early 
developments of conformal field theory (CFT). Starting from a local, classical spin model - e.g. the 
0(n), Potts, Zjv, XY models ~ there are essentially two ways of defining extended interfaces. One 
may consider the polygons in a graphical high-temperature expansion, or the spin interfaces (or 
domain walls) separating regions with different spin values. 

As far as exact solutions are concerned, the situation turns out to be very different for these 
two kinds of interfaces. On one hand, in several models, including the 0(n) and Potts models, the 
high-temperature polygons are described by a Coulomb gas (CG) [Tl[2], i.e. in the scaling limit, they 
become the level lines of a free compact boson CFT. Many geometric properties of these polygons 
can then be computed exactly through their mapping to correlation functions in this CFT. This 
yielded also solid conjectures (and, in some cases, rigorous proofs [3]) about their convergence to a 
Schramm-Loewner Evolution (SLE) model. On the other hand, despite some recent progress [11[5], 
no consistent theory was identified as the scaling limit of spin interfaces, even in well-studied models 
such as the Potts model (except in the particular case of the Ising model), and very few scaling 
properties of these interfaces are actually known analytically. 
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Another current research direction is the study of "extended" SLEs, which are curve models 
determined by two random processes: the driving process of the Loewner chain, and an additional 
process corresponding, in the associated statistical model, to local variables which are not fixed 
locally by the presence of the curve. It was argued [6l [7] that the martingale conditions for this pair 
of processes correspond to null-state equations in an extended CFT. In particular, the presence of 
the additional process affects the relation between the central charge c and the SLE parameter k. 
An important point that needs to be considered in this context, is that in the presence of additional 
variables, the lattice interface no longer satisfies the Markov property, which is an essential ingredient 
of the SLE model. 

The Ashkin- Teller (AT) model is a simple local spin model, with a critical line of constant central 
charge c = 1 and varying exponents. We believe it presents the two features described above: first, 
although it has been much studied with both lattice and CFT techniques, the scaling theory for its 
spin interfaces is not identified, and second, these interfaces leave some fluctuating variables, which 
makes them good candidates for extended SLEs. 

This paper presents two advances in the study of spin interfaces, on the particular example of 
the AT model. First, we introduce some specific boundary conditions (BCs) which ensure that the 
spin interface satisfies the Markov property, even in the presence of an additional variable. We 
then verify numerically that our interface relates to SLE, and for this purpose, we introduce a 
new algorithm allowing the measurement of the left-passage probability [8] in the transfer-matrix 
formalism. Our results also include an accurate numerical determination of the fractal dimension 
df of spin interfaces, also by means of the transfer matrix, which supports a conjecture for df 
stated in [9]. Second, we identify a specific value on the critical line of the AT model, where the 
temperature tends to zero, and at this point, we describe an exact mapping of the AT model onto an 
integrable 19- vertex model. Although this second result does not give direct access to spin interface 
exponents, this represents an alternative way of solving exactly the AT model, and may lead to 
additional exact results on this model. 



2 Spin interfaces in the Ashkin-Teller model 

The Ashkin-Teller model is a two-parameter model with a self-dual, critical line along which the 
central charge is constant c = 1 and the critical exponents vary. Its local correlation functions are 
described by the Z2-orbifold of the compact boson CFT. However, like for the Potts model, this 
CG description only yields results for a specific class of interfaces, which does not include the spin 
DWs. A striking fact illustrating this is the absence of an exact result (except at specific values of 
the couplings) for the fractal dimension of spin DWs along the critical line! 

To be more precise, let us summarize here the current knowledge on spin interfaces in the AT 
model. On the square lattice, the AT model consists of two spin variables aj and tj at each site j, 
taking the values ±1, and with the Boltzmann weight 

where the product is on all pairs of neighbouring sites. The critical line for the square lattice is 
given by the condition: 

sinh2/3 = e-2", (2) 

which we parametrize as 

2 - coth 2/3 = 1 - e"^" = 2 cos irg , with < < 1. (3) 
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To exhibit the Z4 invariance of the model, it is convenient to introduce the complex variables 
Si ; 



l-i 



satisfying 5^ = 1, so that one can write ([T]) as 



W[S] = n exp {f^{S:Sj + S*Si) + I [iS*S,f + iS*S,f] } 



(4) 



e 2 Sj and reflection Sj 



S* 



This Boltzmann weight is invariant under the global rotation Sj , ^ - ^..^ 
A non-branching spin interface in the AT model is specified by the partition into two sets of the 
four possible spin configurations (c, r): 



(+-) 3 : (— ) 



-+)■ 



(5) 



For example, (14|23) denotes the interface between sites with values 1 or 4, and sites with values 2 
or 3. The symmetries of (jH) only leave three possible nonequivalent non-branching spin interfaces: 

• The case (13|24) corresponds to domain walls for the variable r' = err. The associated operator 
has constant conformal dimension h = j throughout the critical line, and hence the fractal 
dimension is |. In a previous work [TO], we constructed a discrete holomorphic parafermion 
■0s of spin s in the AT model, which allowed us to relate the (13|24) interface to SLE(4, p, p), 
where the value of p was derived from the knowledge of s. The problem of the (13|24) interface 
with various boundary conditions was also addressed numerically in [llj and [12], yielding a 
set of conjectured relations to SLE(4, p', p'), and confirming the above results. 

• The case (1|234) is not understood for general couplings, but for g = 3/4 (where the AT model 
coincides with the integrable Fateev-Zamolodchikov Z4 spin model), the associated operator 
satisfies a null-state equation in the Z4-parafermionic CFT, and is argued [13] to have fractal 
dimension j^. 

• The case (12|34) corresponds to domain walls for the variable a. Its fractal dimension is only 
known for a few values of the parameter g: 
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univ. class 


df 


1/4 


0{n = 2) 


3/2 


1/2 


Ising 


11/8 


3/4 


Z4-parafermions 


17/12 


1 


four-state Potts 


3/2 



The fractal dimension cannot be simply obtained from the CG analysis. However, based on 
the above exact values and Monte-Carlo simulations (in the region between Ising and Potts 
points), the following expression was proposed [9] for dj: 



d 



f 



7 g 1 
- + - H 

8 2 85 



(6) 



3 Spin interface with the Markov property 
3.1 Definition 

Let be a domain of the square lattice, and a and b two lattice sites on the boundary. We want to 
find boundary conditions (BCs) such that: 
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• Every spin configuration on il. contains a (12|34) spin interface 7 going from a to b. We denote 
the sites of the path 7 as (70, 71, • • • , 7^), where 70 = a, 7^ = b, and £ is the length of 7. 

• The induced probabihty measure on 7 satisfies the Markov property. 

The Markov property for a probability measure fj.n,a,b on paths means that, for any time s, the 
measure conditioned on (70, . . . , 7^) is the same as the measure defined in 0\(7o, . . . , 7s): 

Vs>0, /^f1,a,6(-|70,---,7s) = /^n\{70,...,7s),7s,6(-)- (7) 

In other words, any piece of curve connected to the boundary should act like the boundary itself, 
when viewed from bulk variables. 

Note that the BCs associated to a (12|34) spin interface - i.e. fixing the boundary spins to be 
in {1, 2} on one side of the boundary between a and b and {3, 4} on the other side - usually do not 
produce a Markov curve as soon as a 7^ 0. Indeed, conditioning on (70, . . . ,7^) correctly fixes the 
a spins on the boundary of Q\{'yo, . . . ,7^), but the resulting measure on spins includes interaction 
factors —OTiTj for neighbouring edges (ij) which cross (70, . . . , 7^), and hence the conditioned mea- 
sure on (7s, . . . , 7^) is not equal to the measure induced by the AT model on r2\(7o, . . . , 7s). We 
thus need to introduce different BCs. 

It is convenient to describe our BCs (see Fig. [T^) on an infinite strip of width L sites 

Sl = {{x,y),x G {1, . . . ,L},y G Z} . 

Spins are fixed to be in {1,2} on the left boundary, and in {3,4} on the right boundary. The 
Boltzmann weights are defined as in ([1]), except that the product is now on neighbouring pairs {ij) 
of the cylinder of width L, i.e. it includes pairs of sites ((l,y), {L,y)). In other words, the BCs are 
such that the a spins live on a strip, and the r spins on a cylinder. This defines a Markov (12|34) 
interface on the strip Sl, since any piece of the interface connected to infinity acts exactly like a 
boundary. 








(a) (b) 

Figure 1: (a) Cylinder of circumference L sites. The thick line is a cut, such that adjacent spins on 
its right have values in {1,2}, and adjacent spins on its left have values in {3,4}. The Boltzmann 
weights include interactions along edges which cross the cut, so that the spins live on the cylinder, 
while the (12|34) interface lives on the strip, (b) The corresponding setting in the infinite plane: 
the domain for the interface is now the plane with a cut {z G C\z < 0}. 
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3.2 Scaling limit 



The above definition may be adapted to any domain 0, formed by the complex plane with an infinite 
cut starting at the origin (see Fig. [lb): the spins are then fixed to {1, 2} on one side of the cut, and 
to {3, 4} on the other side, and we take the AT interaction of the full plane. 

In the scaling limit, the spin interface is expected to become a continuous curve in $7, with a 
conformally invariant distribution. Assuming that the Markov property is preserved in this limit, 
Schramm's principle tells us that the distribution of 7 is SLE^, for some values of k. In the case 
g = 1/2, 7 is an Ising spin interface, and its relation to SLE3 is proved rigorously [3]- 

Note that the above definition is restricted to a particular type of domains, namely those where 
the boundary consists in a cut, i.e. a possibly half-infinite simple curve. This is because we want 
to include interactions across the boundary, which automatically ensures the Markov property. 
Extending this idea to a generic simply-connected domain requires the introduction of a "gluing" 
map / associated to, (f2,a, 6), that is a continuous bijection between the two arcs of the boundary 
defined by a and b. The measure should include an interaction between z and f{z) for all points 
z on the boundary. For instance, in the SLE process, if we denote t the guiding process, ft is a 
one-to-one mapping between {z < at} and {z > at}, which is given by fo{z) = —z at the initial time 
t = 0, and evolves according to the Loewner process. The presence of the process ft certainly affects 
the martingale conditions for correlation functions, which in turn modifies the relation between k 
and the central charge. In particular, in the AT model, this accounts for the fact that the central 
charge remains constant c = 1, whereas k varies as a function of g. 

More formally, one now needs to consider a measure Hn^afi on (7,/), and the Markov property 
reads 

/^n,a,f,(/,-|7[o,t]) = /^n\7rot,,7t,6(/''*V) > where p){z) = \'^^^^ lizedn 

[z itz£-f[o^t]- 

Similarly, under a Loewner process gt in the half-plane, the function ft evolves as 

ft-=9tOp^og^\ (9) 

where /^*^ is defined in ([8]). Assuming that all the relevant information about fluctuating variables 
(in this case, the spins r close to the boundary) is encoded in ft, the Markov property (and the 
associated martingale conditions) should determine a non-standard relation between c and k. 



3.3 Numerical results 

For generic g, we check numerically that the spin interface becomes an SLE in the scaling limit, 
using some known properties of SLE^. First, the fractal dimension of SLE^ is given by 



d 



f 



1 + 



(10) 



Also, the probability that 7 passes to the left of a point (x, y) in the strip Sl is 



1 r(4/K) 
^lefUx) - 2 - ^r(^) 



cot 



/TTX 

It 



2F1 



1 4 3 



n /TTX 



(11) 



In the lattice model, both df and P\cft can be computed numerically by transfer-matrix diagonali- 
sation, yielding two independent determinations of k. The fractal dimension is obtained from the 
two- leg watermelon conformal dimension X2 as dj = 2 — X2, where X2 is given by the dominant 
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eigenvalue of the transfer matrix on the cyhnder in the sector where two interfaces are forced to 
propagate along the cylinder. Moreover, we observe a better convergence when the model lives on 
the triangular lattice, so we chose this setting for the determination oi dj. The probability Pieft is 
computed from the components of the dominant eigenvector of the transfer matrix on the strip, as 
described in the Appendix. 

The results for dj are shown in Fig. [2l and are in very good agreement with the conjecture ([6]), 
on a large portion of the critical line. Close to g = 1, marginal operators typically arise, and cause 
very strong finite-size effects, which explains the poorer convergence in this region. The validity 
of dni) in this region was already checked by Monte-Carlo simulations in [9]. 

For the computation of Picft, we use the value g = 1/3, where the associated model is simpler, and 
hence the space of states on which the transfer matrix acts has smaller dimension. The numerical 
results on Pieft, and the corresponding estimate of k, are shown in Fig. [3] and Tabled! They match 
closely with the results on dj, which is a strong indication that the curve distribution is SLE^^xq/s- 




Figure 2: (Color online) Numerical determination of the fractal dimension df for the (12|34) interface 
in the triangular-lattice AT model. The dotted curve represents the conjectured value ([6]). 



4 The zero-temperature point 

The value g = 1/3 on the self-dual line ^ corresponds to the limit 
The Boltzmann weights then take the form 



(12) 
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Figure 3: (Color online) Numerical determination of the left-passage probability Picft in the AT 
model for g = 1/3. The full line represents Schramm's formula, with the value k = 10/3 from the 
conjecture dH). 



where M is the number of edges in 17. The corresponding a and r domain walls thus form a two- 
colour loop model on the dual lattice, where a black (resp. gray) loop represents a a (resp. r) 
domain wall. Each edge of the dual lattice bears at most one colour of loop, and the weight of a 
loop configuration is K^, where K = 1/2, and I is the total length of the loops. For completeness, 
we depict in Fig. H] all the possible loop configurations around a vertex. 
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Figure 4: Vertices of the two-colour loop model. 



Let us now expose an exact equivalence of this two-colour loop model of AT spin interfaces with 
an integrable 19-vertex (19V) model. As a first stage, we map the 19V model onto a dilute Brauer 
model with loop weight n = 2. The vertices of the latter model are shown in Fig. [5j Each of the 
Brauer loops can be oriented clockwise or anti-clockwise independently: this leads to a 19V model 
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Table 1: Numerical estimates of n from the left-passage probability Pieft for ^ = 1/3 (see data in 
Fig. ED. 
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where we have used the standard notations (wi, . . . ,1^19) as in 
weights [H] of the 19V, we obtain the dilute Brauer weights: 



t = sin//sin2/i — siniisin(/i — li) 

Ul = sin 2/i sin(/x — n) 

U2 = sin 2/i sin u 

V = sin u sin(/x — u) 

Wl = sin 2/i sin(/x — n) cos n 

11)2 = sin2/isinncos(^ — n) 

X = — cos 2/i sin n sin(/x — n) . 



(13) 
(14) 
(15) 
(16) 
(17) 
(18) 
(19) 

If we take the FZ integrable 



(20) 



It turns out that setting = ^ and u 
appear in even number) yields 



fi/2, and letting v — t- —v (because v configurations always 
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(21) 



The dilute Brauer model is, in turn, mapped to our two-colour loop model in the following way. 
Each of the Brauer loops can be coloured black or gray independently. Then, consider a vertex 
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whose four adjacent edges bear the same colour. Since the coloured loops have weight unity, the 
overall Boltzmann weight of the configuration is the same for vertices wi,W2 and x. Thus the dilute 
Brauer model is equivalent to the two-colour loop model provided the following relations hold 

ai 

02 = • • • = as 
ag = • • • = ag 

o-io = • • • = ai3 
ai4 = ai5 

am + 020 = wi9 + W23 
an = ais 
021 = a22 

The integrable weights (12ip satisfy these relations, and hence we have proved the equivalence. 

For the parameters = ^ and u = /i/2, the 19V model is described [16] in the continuum limit 
by a c = 1 compact boson CFT with the spectrum of conformal dimensions 

Xn,m = i^n^ + 7T-"^^ ■, n,mel? . (23) 
27r 2/x 

It turns out that the excitations described by this spectrum describes are not directly related to 
those of the colored loops. In particular, they do not include the conjectured two-leg exponent 
X2 = 7/12. However, the mapping to a 19V model may be exploited to obtain new exact results on 
the AT model at zero temperature. 



t 

Ul 

U2 
V 

zizX 

wi + -0:2 + X 

Wl 
W2 



(22) 



5 Conclusion 

In the present work, we have presented a set of new results on spin interfaces in the AT model. 

By a numerical study based on transfer-matrix diagonalisation, we have given a precise numer- 
ical determination of the fractal dimension dj, in a parameter regime that could not be accessed 
previously by a Monte-Carlo simulation [9]. Our results on df bring more support to the conjec- 
tural expression ^ introduced in [9]. Moreover, we have shown, for the first time, how to adapt the 
transfer-matrix method to measure numerically the left-passage probability Pieft, which is a natural 
quantity in the SLE literature. 

On the analytical side, we introduce a general way of choosing BCs which preserve the discrete 
Markov property of a path, even in the presence of additional degrees of freedom which are not 
fixed locally by the path. This construction is quite general, and may be extended to any lattice 
model with interfaces of this sort. The left-passage probability Picft is then computed numerically 
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with these "Markovian" BCs, and the results for Picft confirm that the resulting measure on curves 
becomes SLE in the continuum limit. 

We have also identified a specific point on the critical line of the AT model, the point where all 
the couplings diverge, where there is a mapping of the AT model onto an integrable 19- vertex model. 
At this point the interfaces of the a and r spins live in different bonds and so the corresponding 
two-color loop model has a simpler form. At this particular point one can also find a dilute Brauer 
model representation of the critical spin interfaces of the AT model. Since the loop model has 
a simpler form at this particular point, we believe it is a good place to look for a Coulomb gas 
representation of the AT spin interfaces, and extend it to the critical line afterwards. 

Our work presents conclusive numerical arguments that the spin interfaces in AT model are 
conformally invariant and their fractal dimension follow the formula ([6]), motivating further the 
search for an analytic argument which would account for ([6]). This is a simple example of the 
more general notion of spin interfaces, for which very few analytical results are available. On top 
of the recent methods for the study of lattice models, like discrete holomorphic parafermions, our 
construction of "Markovian" BCs seems an interesting starting point for this search. 



Appendix: transfer-matrix algorithm for the left-passage probabil- 
ity 

In this Appendix, we present a numerical procedure for the computation of the left-passage proba- 
bility in a loop model defined on the infinite strip of width L sites. Let us describe how our method 
works on a simple loop model, such as the square-lattice 0(n) loop model (defined by the vertices 
in Fig. [51 with x = 0). The generalization to other loop models, such as the two-colour loop model 
studied in this paper, is straightforward. 
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Figure 6: a) Some basis elements for L = 5. b) An example of the gluing operation involved in the 
definition of (•,•). 



The basis elements a generating the space of states Vl are labelled by planar pairings of the 
points {1, . . . , L}, with one point paired with as depicted in Fig. [6l The space Vl is also equipped 
with a scalar product (•,•), defined as follows. We say that two pairings a,/3 are compatible if 
they have the same subset of paired points, and we denote a/3 the figure obtained by gluing the 
TT-rotation of a to /3. Then 

^ _ in*"^"^'"^ i°°P^ '"^ if Q and /3 are compatible, ^^^^ 
1 otherwise. 
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Let us recall briefly how the transfer matrix T acts on Vl |15] . We take the convention that T acts 
from bottom to top, and adds one row to the system. When starting from a pairing configuration 
1/3), the addition of one row gives rise to a linear combination of pairings X^Q,Ta/3|a), where the 
matrix element T^/s is given by the product of Boltzmann weights inside the row, including a factor 
n for each loop that got closed in the process. 

We denote the dominant eigenstate of T in the sector with a propagating path as 

|V)=^^a|a). (25) 
a 

If we cut the system into two half-infinite strips above and below a given horizontal section, then 
the L dangling edges of the top half have a pairing a, and those of the bottom half have a pairing 
/3. If we assume that the local Boltzmann weights in T are invariant under a rotation of angle vr, 
then the probability of a pair (a, f3) is 

F(a,/3) = ^"^^j"'^\ where Z = Y,Mp{a,(3) (26) 

a,l3 

For any pair (q;,/3), we denote ^a,i3 the open path appearing in a/3. For a given position 
j G {0, . . . , L}, we introduce the modified scalar product (•, 

I n\ I if a and /3 are compatible and 7^ r passes between and j, . 

{a,l3)j:={ . (27) 

I (J otherwise. 

This allows us to express the left-passage probability Pioft as 

^ieft(j) = ^^^V'a^ (a,/3)i. (28) 

Our method consists in the following steps 

1. Determine the eigenstate ()25p by the power method (iteration of T on a random initial vector). 

2. Compute (j28|) by performing the double sum on configurations. 

Since the computation of {a,/3)j takes 0(L) operation, the overall complexity of step 2 is 
0[L(dim Vl)^]- The cost in memory is 0(dim Vl), since only a finite number of vectors are stored 
in memory, as usual in transfer-matrix algorithms. 
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